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Abstract 

Using the technique of quantum transport of the interfering pair we examine the Hanbury- 
Brown-Twiss (HBT) interferometry signatures for the particle-emitting sources of pions and kaons 
produced in the heavy ion collisions at GSI-FAIR energies. The evolution of the sources is described 
by relativistic hydrodynamics with the system equation of state of the first-order phase transition 
from quark-gluon plasma (QGP) to hadronic matter. We use quantum probability amplitudes in 
a path-integral formalism to calculate the two-particle correlation functions, where the effects of 
particle decay and multiple scattering are taken into consideration. We find that the HBT radii 
of kaons are smaller than those of pions for the same initial conditions. Both the HBT radii of 
pions and kaons increase with the system initial energy density. The HBT lifetimes of the pion and 
kaon sources are sensitive to the initial energy density. They are significantly prolonged when the 
initial energy density is tuned to the phase boundary between the QGP and mixed phase. This 
prolongations of the HBT lifetimes of pions and kaons may likely be observed in the heavy ion 
collisions with an incident energy in the GSI-FAIR energy range. 
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I. INTRODUCTION 



One of the central goals of high enery heavy ion collisions is to find and quantify the QCD 
phase transition from hadronic matter to quark-gluon plasma (QGP). Under the assumption 
of first-order phase transition there is a mixed phase of the QGP and hadronic gas. In the 
absence of pressure gradient, a slow-burning fireball is expected when the initial system is at 
rest in the mixed phase, and this leads to a considerable time-delay of the system evolution 
. It is therefore of interest to investigate the time-delay signatures for the 
first-order phase transition. 

Two-particle Hanbury-Brown-Twiss (HBT) interferometry is a useful tool for detecting 
;he space-time structure of particle-emitting sources in high energy heavy ion collisions 



3, y, y, 



101 ] . For the first-order phase transition between the QGP and hadronic matter, 



the time-delay of the system evolution may prolong the emission duration of particles and 
lead to unusually large HBT lifetime^as compared to a crossover transition or a hadron gas 
without the QCD phase transition . It is known that the phase transitions 

occurring in the heavy ion collisions at RHIC and top SPS energies are crossover and in small 
baryon density regions. Moreover, at AGS energies the systems are almost in hadronic phase 
although with higher baryon densities. The future Facility for Antiproton and Ion Research 



(FAIR) at GSI with heavy ion beams from 2 - 45A GeV will provide an 
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explore the first-order QCD phase transition at high baryon densities jll . 
Investigating how the HBT results variate with the incident energy and what the HBT 
signatures the first-order phase transition is possessed of in the GSI-FAIR energy range is 
thereby the subject of this work. 

In high energy heavy ion collisions the final particles include the contributions of direct 
production and excited-state particle decay. Also, the particles will subject to the multiple 
scattering with other particles when they propagate in the system up to the thermal freeze- 
out. In Ref. a HBT analysis technique with quantum transport of the interfering pair 
was developed to investigate the effects of particle decay and multiple scattering on the 
extracted HBT radii in the heavy ion collisions at AGS and RHIC energies. This HBT 
analysis technique allows one to follow the trajectories of the test particles after emission up 
to the thermal freeze-out in model calculations. It is suitable for examining the space-time 
geometry of evolving sources in detail. In this study we will use relativistic hydrodynamics 
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with an equation of state of first-order phase transition to describe the evolution of the 
particle-emitting sources produced at GSI-FAIR energies. We will use the HBT analysis 
technique with quantum transport of the interfering pair to examine the HBT radii and 
lifetimes of the sources for different initial energy densities. Because the central heavy ion 
collisions at GSI-FAIR energies are almost full stopped, we assume that the systems are at 
rest initially in the center-of-mass frame and with spherical shape for simplicity. 

As compared with pion HBT interferometry kaon HBT interferometry may present more 
clearly the source space-time geometry at emission configuration, because kaons can escape 
easily from the system after hadronization and therefore seldom be affected by the multiple 
scattering and particle decay. By comparing the results extracted from the two-pion and 
two-kaon HBT interferometry, we find that the multiple scattering and particle decay lead to 
a larger radius of the pion source than that of kaon's. However, both the HBT lifetimes of the 
pion and kaon sources have significant increase when the initial energy density approaches 
to the boundary between the QGP and mixed phase. 

The paper is organized as follows. In section II we describe the model of equation of 
state (EOS) used in our calculations. The adiabatic paths for the EOS and the space-time 
configuration of the evolving system are also discussed in this section. In section III we 
give a brief description to the quantum probability amplitudes in a path-integral formalism 
of HBT and the two-pion and two-kaon HBT results for the evolving sources produced in 
the collisions at GSI-FAIR energies. Finally, the summary and discussion are presented in 
section IV. 



II. EQUATION OF STATE AND SYSTEM EVOLUTION 

For non-dissipative ideal fluid hydrodynamics is defined by the local conservations of 
energy-momentum and other conserved quantities (e.g. entropy, baryon number, and 
strangeness) Q- To solve the hydrodynamical equations one needs the EOS whi ch g ives 
the relation among the thermodynamical quantities in the conserved equations 18|, ll9|. In 
what follows we discuss our EOS model and system evolution described by hydrodynamics 
and the EOS. 
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A. EOS model 



In our EOS model the QGP phase is described by an ideal gas of gluons, u, d, s quarks and 



20 



2l| 



antiquarks, with the constant vacuum energy B associated with QCD confinement 
The pressure, energy density, and the conserved charge density in the QGP phase are given 
by 



P 



Q 



^Pi(T,Hi) - B 



e Q = J2 £ i( T ^i) + B 

i 

n Q A = rii(T, fj,i) , 



(1) 
(2) 
(3) 



where Pi(T, fii), £j(T, /ij), and nj(T, /Zj) are the pressure, energy density, and number density 
of particle species i in the ideal gas with temperature T and chemical potential {/ij}, A^ is 
the conserved charge number of the particle species i. In our calculations we use the quark 
masses m u = = 5 MeV, m s = 150 MeV and the bag constant B = (235 MeV) 4 21]. 



For the hadronic phase we adopt the excluded volume model 
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231 ] and consider 



the particles ir, K, N, A, E, A, and their antiparticles. The pressure, energy density, and 
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22j,|23j 



(4) 

(5) 
(6) 



the conserved charge density in the hadronic phase are given by 

i 

where 

fn= fXi-V Q p H , (7) 

Vo = (l/2)(47r/3)(2a) 3 is the excluded volume which is assumed to be the same for all 
hadrons with a = 0.5 fm 21]. 

For the first-order phase transition, there are Gibbs relationships in the mixed phase of 
the QGP and hadron gas. We have T Q = T H , /ijv,A — 3/i u , //a,s — 2/i u + fi^+^o^- = 0, 
^k+,k° = Hu — Us, • and 

P M = P Q (T, [in, Us) = P H (T, ft u , »s) , (8) 



a 



e Q (T, fi u , n s ) + (1 - a) e H (T, /i u , /z s ) 



(9) 



nf = a nj(T, fi u , fi 8 ) + (1 - a) n H A (T, fi u , fi 8 ) , (10) 

where \x u and /x s are the chemical potential of u and s quarks, and a = Vq/V is the fraction 
of the volume occupied by the plasma phase. The boundaries of the coexistence region are 
found by putting a = (the hadron phase boundary) and a = 1 (the plasma boundary). 

Using the thermodynamical relations of ideal gas one can get other thermodynamical 
quantities, such as entropy density s, in the QGP, hadronic, and mixed phases from Eqs. 
~ ©], [fllj) - ©], and [(J8l)- (jl0|) ]. and get numerically the EOS for solving the hydrodynamical 
equations. 



B. Adiabatic paths 

In our model calculations the system evolves from a thermalized initial state to final 
freeze-out. In the absence of dissipation, the entropy of the system is conserved during 
evolution. On the other hand, the baryon number is also conserved. So the ratio of their 
local densities %/s is a constant. For our EOS model we show in Fig. 1 the adiabatic cooling 
paths of the systems with ns/s =0.08 and 0.06, which correspond to the incident energies 



about 10 and 30 AGeV, respectively 
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24j . The dotted line in Fig. [T] is the transition 



curve between the QGP and hadron gas. The mixed phase is on the transition curve from 
the endpoint of the QGP branch (point 1 or 1') up to the beginning of the hadronic branch 
(point 2 or 2'). The non-trivial zigzag shape of the trajectories indicates that the system 



has a re- heating in the mixed phase 23j, |25[ . The reason is that at a certain point (T, /x) on 
the phase-transition curve, the number of degrees of freedom and hence the specific entropy 
in the plasma phase are larger than the corresponding values in the hadronic phase. Hence 
the temperature must increase during hadronization to conserve both the total entropy and 
baryon number simultaneously [3]. 

In Fig. [2] we show the hydrodynamical relevant relation, p/e, for hb/s =0.06 and 0.08. 
At the boundaries between the hadronic and mixed phases (2 and 2') the ratio p/e has 
maximums. The minimums of the ratio, named the "softest points" , are at the boundaries 
between the mixed phase and QGP (1 and 1'), corresponding to e = e MCi =1.83 and 1.90 
GeV/fm 3 for %/s = 0.06 and 0.08, respectively. 
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FIG. 1: Adiabatic paths for tib/s = 0.06 and ub/s = 0.08. 
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FIG. 2: The ratio of pressure to energy density p/e for hb/s = 0.06 and ub/s = 0.08. 
C. System evolution 

equations for 



After knowing system EOS we can obtain the solutions of hydro dynamica 

by using the HLLE scheme 



the certain initial conditions 
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30| and Sod's 



operator splitting method 31] . Because the central heavy ion collisions at GSI-FAIR energies 
are almost full stopped, we assume that the initial system is at rest in a sphere with a constant 
energy density e°. For ub/s — 0.06, the incident energy is about 30 AGeV 



23 



|2J. We 

investigate the system evolution with the initial energy densities e° = 4.12 GeV/fm 3 > e MCi 
1.83 GeV/fm 3 . The corresponding initial temperatures are T° =180 and 142 



and e c 



-MQ 



MeV. Meantime, the corresponding initial baryon chemical potentials are fi° B = 3/i° =990 
and 780 MeV. For the system with ns/s = 0.08, the corresponding incident energy about 10 
AGeV enable only the initial energy density to approach the region of the QGP boundary 
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FIG. 3: The space-time contours of energy density for the systems with ng/s = 0.06 (a,b) and 
ub/s = 0.08 (c,d), for the initial energy density e° > e MCl , e° = e MCl , and e° = e HM . 
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241 ] . In this case we calculate the system evolution with the initial energy densities 
e° = e MQ = 1.90 GeV/fm 3 and e° = e UM = 172 MeV/fm 3 , where e UM is the energy 
density at the boundary between the hadronic and mixed phases. The corresponding initial 
temperatures are 132 and 152 MeV. The corresponding initial baryon chemical potential are 
960 and 480 MeV. 

Figure [3](a) and (b) show the space-time contours of the local energy densities at e MCi , 
£ HM , and e TFO for e° > e MQ and e° = e MCi for the system with Ub/s = 0.06, respectively. 
Here r and e TFO are the initial system radius and the energy density at the thermal freeze- 
out. One can see that the system takes more time evolving through the mixed phase (from 
£ mq £ Q £ hm-j ^jjgjj that through the pure QGP phase (from e° to e MCi ) or the pure hadronic 
phase (from e RM to e TFO ). The duration of the evolution through the mixed phase is larger 
when the initial energy density is at the soft point. Figure M^c) and (d) show the contours 
at e HM and e TFO for e° = e MCl and e° = e HM for the system with ns/s = 0.08, respectively. 
For e° = e M Q, because the soft point for ns/s = 0.08 is higher than that for tib/s = 0.06 
(see Fig. [2]), the evolving time through the mixed phase for ub/s = 0.08 is short than that 
for ns/s = 0.06. For e° = £ HM there is only hadronic phase. The evolution is fast because 
of the larger p/e in the hadronic gas. In our calculations the energy density at the thermal 
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freeze-out is taken to be 45 MeV/fm 3 32j, which corresponds to the thermal freeze-out 
temperatures 110 and 100 MeV for %/s =0.06 and 0.08. The initial system radii are taken 
to be 6 fm. 



III. HBT INTERFEROMETRY WITH QUANTUM TRANSPORT OF THE IN- 
TERFERING PAIR 

The two-particle Bose-Einstein correlation function C(ki, k 2 ) is defined as the ratio of 
the two-particle momentum distribution P(ki,k 2 ) to the the product of the single-particle 
momentum distribution P(ki)P(k 2 ). For an evolving source, using quantum probability 
amp litudes in a path-integral formalism, P(ki) (i = 1,2) and P(ki,k2) can be expressed as 



13, 



26 
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34 



35] 



P{ki) = j d i xp{x)e- 2Jm ^ x) A 2 {Kx) , (11) 

P(k u k 2 ) = J d 4 Xl d 4 x 2 e -2W.(*i) e -2Zni*.0w) 

xp(x 1 )p{x 2 )\<!>(x 1 X2;k 1 k 2 )\ 2 , (12) 

where p(x) is the four-dimension density of the particle-emitting source, A(nx) is the mag- 
nitude of the amplitude for producing a particle with momentum k at x, e - 2Xm< t>s(. x ) [ s the 
absorption factor due to multiple scattering, and Q(xix 2 ; k\k 2 ) is the wave function for the 
two identical boses produced at X\ and x 2 with momenta /ti and k 2 , and detected at x&_ or 
Xd2 with momenta k\ and k 2 , respectively. 

In our HBT calculations final identical kaons are considered to be emitted thermally 
from the space-time hypersurface at e HM and keep freeze-out after emission. However, final 
identical pions (for example tt + ) include the primary pions emitted from the hypersuface 
at e HM and the secondary pions from the "excited-state" particle decays during the system 
evolving in hadronic phase until to the thermal freeze-out. The four-dimension density of 
the pion source can be expressed as 17] 



p(x) = n n (x)8(t - r°) + ^2 Dj^ w nj(x) , (13) 

where rii(x) and r° are the particle number density and the hadronization time in local 
frame, Dj^ n is the product of the decay rate in time and the fraction of the decay dj^. 
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For example, 



Ta x | and D^o^o 



7r + 7r ) x 1, where v r is 



the relative velocity of the two colliding pions and the cross section <7(7r°7r° 



7T 1 7T 



IS 



7T 7r° 



17J. In calculations we neglect the 



equal to the absorption cross section of n + n~ 
contributions of A and S decays to the final pions, because most of them exist until the 
thermal freeze-out. 

When a test pion propagating in the source it will subject to multiple scattering with 



the medium particles in the source. Based on Glauber multiple scattering theory 



36|], the 



absorption factor due to the multiple scattering in Eqs. (TTT|) and (TT2"1) can be written as 



[13, 



26 



33, 



34 



35] 



-21m 4> s (x) 



exp 



^2 °"abs(vri) ni(x')jd£(x' 



(14) 



where Y2i means the summation for all medium particles except for the test pion along the 
propagating path d£(x'), a^^Tii) is the absorption cross section of the test pion with the 
particle species i in the medium, and Xf is the freeze-out coordinate. In calculations we only 



consider the dominant absorption processes for the identica 
of 7T + 7r~ — > 7r°7r° and it + N — > A for 7r + , as we did in Ref. 



pions, for example the reactions 



17|. 



In HBT analyses the variables usually used are the Pratt-Bertsch variables, g out , g S ide, and 



ong 



23, 



38(. Here g out and g s ;d e are the components of the relative momentum of identical 



particle pair in the directions parallel and perpendicular to the total transverse momentum 
of the pair, and gi ong is the relative momentum along the beam direction of collision. Using 
the Pratt-Bertsch variables the HBT radius in the side-direction, -R S id e , reflects the transverse 
size of the source. However, the HBT radius in the out-direction, R nui , is related to not only 
the source size, but also the source expanding velocity and lifetime 37], |38(. So a detailed 



joint analysis of R out and -R S ide as a function of transverse momentum of the pair may also 



provides the information of source dynamics 



37 



38j |. Motivated by investigating the 
source lifetime directly and clearly, we use the variables q — |ki — k 2 | and go = \E\ — E 2 \ 
and the simple Gaussian fitting formula 



C(g,g ) = l + Ae-o 2i?2 -^ 2 



(15) 



where R, r, and A are the source HBT radius, lifetime, and chaotic parameter. 

From Eqs. (TTTD (fT4|) ^ 
function for each (g, q ) bin 



From Eqs. (FlT!) - (fl4j) we can construct numerically the two-particle HBT correlation 

39I I40I . Figure 0] shows the two-particle correlation 
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FIG. 4: The two-particle correlation functions of kaons and pions for the evolving systems as shown 
in Fig. El 

functions C(q, go <15MeV/c) for the evolving sources with %/s =0.06 and 0.08, and for 
e° > e M Q, e° = e M Q, and e° = e nM . The symbols a and • are the two-kaon and two-pion 
correlation function results. For comparison, the symbols □ present the two-pion correlation 
functions calculated with the pions emitted from the thermal freeze-out (TFO) configuration. 
Table [I] gives our HBT fitted results. One can see that for the same initial conditions, the 
two-kaon HBT radius is smaller than that of the two-pion's. It is because that the kaons are 
emitted earlier and from smaller system configuration. The multiple scattering and particle 
decays during the source evolving in hadronic phase increase the HBT radii of the pion 
sources. Because the pion sources have largest configuration at TFO, the corresponding 
HBT radius is the largest for the certain initial energy density. 

For the system with %/s = 0.06 we can see that for the initial energy density e° > e MCi 
the two-pion HBT lifetimes are larger than that of two-kaon's. However, when e° drops 
to the soft point e M( ^ the HBT lifetimes for the pions and kaons increase significantly and 
almost are the same, while the corresponding HBT radii decrease. The reasons are related to 
their different source space-time geometries [see Fig. [3] (a) and (b)] and expanding velocities 
for the two kinds of initial conditions. When the initial energy density is at the soft point, 
the source has the smallest expansion and correspondingly the smallest average spatial size 
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TABLE I: The HBT fitted results. 



2K 2ir 2tt(TF0) 



(a) R = 3.25 ±0.12 fm R = 4.90 ± 0.20 fm R = 6.95 ± 0.30 fm 
e° > e MQ r = 5.35 ± 0.36 fm/c r = 7.31 ± 0.54 fm/c r = 7.84 ± 0.74 fm/c 
ni3/s = 0.06 A = 0.81 ±0.05 A = 0.83 ± 0.05 A = 0.90 ± 0.06 

(b) R = 2.61 ± 0.09 fm R = 3.70 ± 0.14 fm R = 5.20 ± 0.19 fm 
£ o = £ mq r = 12 50 ± 68 fm/ / c T = 12 53 ± 58 fm / c r = 12 43 ± o.65 f m / c 

n B /s = 0.06 A = 0.91 ±0.05 A = 1.04 ± 0.05 A = 1.03 ±0.05 

(c) R = 2.35 ± 0.07 fm R = 3.70 ± 0.10 fm R = 4.98 ± 0.16 fm 
£ o = £ mq> r = 9.58 ± 0.36 fm/c r = 9.52 ± 0.31 fm/c r = 9.13 ± 0.35 fm/c 
n B /s = 0.08 A = 0.98 ±0.05 A = 1.04 ±0.04 A = 0.95 ± 0.04 

(d) R = 1.66 ± 0.03 fm R = 2.38 ± 0.06 fm R = 3.54 ± 0.04 fm 
e° = e HM r = 2.33 ± 0.08 fm/c r = 2.43 ± 0.10 fm/c r = 1.73 ± 0.14 fm/c 
n B /s = 0.08 A = 0.93 ±0.01 A = 1.02 ±0.02 A = 1.00 ±0.01 



and largest evolution time from the initial state to the hadronization. Because there is no 
the influence of the system evolution on kaons after hadronization, the increase of the HBT 
lifetime of the kaons reflects the prolongation of the system evolution in the mixed phase. 
In Table I, the results of chaotic parameter A for e° > e MCi are obviously smaller than unit. 
This is mainly because that the particle-emitting sources in this case are much different from 
Gaussian distribution. 

For the system with ub/s = 0.08, the two-pion and two-kaon HBT lifetimes for e° = e MCi 
are much larger than those for e° = e RM . Also, the HBT radii for e° = e MCi are larger 
than the corresponding results for The main reason for these results is that the 

space-time configuration for e° = e RM is much smaller than that for e° = e MCl [see Fig. 
[3] (c) and (d)]. Our HBT investigations indicate that both for %/s =0.06 and 0.08, the 
maximums of the HBT lifetimes appear when the initial energy density reaching at the soft 
points. Because the ratio p/e at the soft point e MCl for ub/s = 0.08 is higher than that for 
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ub/s = 0.06 (see Fig. [2]), the maximal HBT lifetimes for ns/s = 0.08 are smaller than those 
corresponding lifetimes for ns/s = 0.06. 

For the systems with the Ub/s values between 0.06 and 0.08, we also find that their HBT 
lifetimes increase significantly when e° approaches at the corresponding soft points (between 
1 and 1' in figure [2]). The maximums of the HBT lifetime are about 10 fm/c much larger 
than the results for e° > e M( ^ and e° = £ HI 1_L in Table I (a) and (d)]. Based on the evolving 



trajectories calculated by hydrodynamics 
staying in the mixed phase of the QGP and hadronic gas wi 
energies are between 10 and 30A GeV (see the figure 19 of [24 



241 ]. the events with the initial thermalized states 
1 happen when the incident 
j, the beginnings of the bold 



parts of the trajectories for the 10 and 30v4 GeV are at the two sides of the transition region, 
respectively). Correspondingly, the ratio of %/s is about between 0.08 and 0.06 (see the 



figure 18 of 



241] ). From our model calculations, the maximums of the two-pion and two-kaon 



HBT lifetimes will be observed simultaneously in GSI-FAIR energy range when the initial 
energy density is tuned to the soft point. 

IV. SUMMARY AND DISCUSSION 



Recently the heavy ion collisions at the energies between AGS and the top-energy SPS 



attract special attention, for example the SPS and RHIC low energy programs 
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42 



43 



and the project of GSI-FAIR 
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15 



44] 



In this energy range it is expected that 



the heavy ion collisions may produce the QGP with high baryon density, which is differen 
from that have been observed in RHIC and top-energy SPS experiments. Correspondingly, 
the phase transition from the high-baryon-density QGP to hadronic matter is the first-order 
transition, which will lead to a system evolution much differen from the crossover transition 
happened in the low baryon density region at RHIC and top SPS energies. 

Using the technique of quantum transport of the interfering pair we examine the two- 
pion and two-kaon HBT interferometry for the particle-emitting sources produced in the 
heavy ion collisions at GSI-FAIR energies. We use relativistic hydrodynamics with the EOS 
of first-order phase transition between the QGP and hadronic gas to describe the system 
evolution. The two-particle HBT correlation functions are calculated with the quantum 
probability amplitudes in a path-integral formalism, where the effects of particle decay and 
multiple scattering are taken into consideration. We find that both the HBT radii of pions 
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and kaons increase with the system initial energy density. The particle decay and multiple 
scattering lead to the larger HBT radii of pions than the corresponding HBT radii of kaons. 
The HBT lifetimes of the pion and kaon sources are sensitive to the initial energy density. 
They are significantly prolonged when the initial energy density is tuned to the soft point. 
Our model calculations indicate that this significant prolongation of the HBT lifetimes of 
pions and kaons will be observed in the heavy ion collisions at GSI-FAIR energies. 

As a useful space-time probe HBT interferometry has been extensively used in high energy 
heavy ion collisions. However, there are still some open problems on HBT analysis technique 
and the understanding of HBT results. The HBT measurements at RHIC indicate that the 
values of the ratio of the transverse HBT radii R out to -R S ide are smaller than those from the 



hydrodynamical calculations 



45 



46 



47], 



forth to explain the RHIC HBT puzzle 



21, 



481 ]. Various models and techniques have been put 
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52 



53 



54 



55 



56 



57 



58, 



59, 6o], 



6l\ . At GSI-FAIR energies, the heavy ion collisions are almost full stopped. We used an 
approximation of spherical evolving sources and assumed that the initial states are static 
and uniform in our calculations. It would be interesting to consider more reasonable evolving 
sources and study the effect of initial conditions on the HBT results in future investigations. 
Also, a systematical investigation to HBT interferometry in different energy ranges will be 
of great interest. 
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